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Introduction 



The study of the effect of pressure on materials is fundamental to a range of prob- 
lems in physics (e.g., Ref. [1]). With increasing pressure, the distances between atoms 
normally decrease, which leads to the increase of interparticle interactions. Thus, by 
changing the pressure, one can tune the physical properties of materials and observe cor- 
responding changes in the vibrational, electronic and magnetic excitations. The gradual 
changes of interatomic distances and interaction in some cases yield to radical transfor- 
mations that are unique to high pressure. The behavior of compressed materials can 
be studied under static conditions to pressures above 300 GPa (or 3 megabars). The 
associated large decrease in volume provides a wealth of opportunities to observe basic 
changes in bonding character, magnetic or electronic structure or chemical state. Such 
experiments have become possible with the development of ultrahigh-pressure diamond- 
anvil cell techniques (e.g., Ref. [2]). In particular, high-pressure spectroscopy provides 
crucial and often unique information on bonding properties and excitations of metals, 
semiconductors and superconductors [3, 4]. Ultrahigh-pressure research, by its very na- 
ture, imposes substantial requirements on experimental technique. The extremely small 
amount of the material (e.g., down to picoliter volumes) brought to megabar pressures 
demands extremely sensitive and versatile equipment. Thus, the quality of information 
obtained critically depends on the development of new techniques and improvement of 
the existing methods. 

Raman studies of metals is a challenging task because only a very thin surface layer 
of the sample (skin-depth of <1000A) is typically probed by the exciting laser radiation. 
This makes the requirements for the experimental technique even more severe. Raman 
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studies in metals were began in 1966, with the pioneering work of Feldman et al. [5]. 
They were able to detect very weak signals from phonon modes in Be and AuAl 2 . The 
first work on metals under pressure were published in 1992 by Olijnyk [6]. The high- 
pressure phases of Si and Ge were studied to 50 GPa by Raman measurements of the 
phonon modes, and the results were in reasonable agreement with theoretical calculations 
[7, 8, 9]. Such studies have now become more feasible, reliable, and definitive with the 
development of new CCD array detectors and optical systems [10, 11], because they give 
tremendous advantage in comparison to single-channel and diode-array techniques. 

Raman spectroscopy of metals has been of growing interest, in part because it is 
a very informative way of characterizing high-temperature superconductors (including 
determination of the superconducting gap; e.g., Refs. [12, 13]). The possibility of using 
Raman spectroscopy under pressure for studies of these materials has been demonstrated 
[14]. Several applications of the technique to metals have been reported [15, 16, 18], but 
its more extensive use has been limited in pressure range and by the need to detect very 
weak signals. We applied new holographic transmission optics Raman techniques [19] for 
high-pressure studies to obtain substantially improved performance, which allowed us to 
extend Raman measurements in metals to the megabar pressure range [20]. 

2. Experimental 

2T. Raman technique. - The recent improvements in Raman spectroscopy techniques 
stemmed from the development of extremely sensitive array detectors (CCD) and fast 
imaging spectrometers [4, 11] coupled with and incorporating holographic transmission 
optics. The latter has already found a wide application as notch filters, which efficiently 
block undesirable Rayleigh-scattered light in Raman systems used for routine measure- 
ments [21]. This obviates the need for cumbersome low-transmission double subtractive 
monochromator filter systems for measurements above 50 cm -1 and increases substan- 
tially the throughput of the system. The obvious drawback of such systems in its com- 
mercial realization is the need of a separate set of optics for each excitation wavelength 
(e.g., lines of an Ar-ion laser). 

In the case of high-pressure spectroscopy, flexible changes in both the excitation wave- 
length and the spectral range are essential because of the need to choose the optimum 
excitation for both ruby fluorescence and Raman measurements. This requires conve- 
nient tuning of the central wavelength of the laser bandpass filter and the position of the 
notch of the laser rejection filter. We have developed a simple universal design for Ra- 
man/fluorescence measurements with capabilities of collecting low- and high- resolution 
spectra and exciting spectra in different geometries using a variety of laser wavelengths, 
which can be easily selected by tuning holographic optics (fig. 1). The spectral and spa- 
tial purity of the laser beam is greatly improved by use of a bandpass filter (Kaiser Optical 
Systems) and a spatial filter, which is also used as a beam expander. The collimated 
laser beam is focused to the sample in either backscattering geometry with a neutral or 
dispersive beamsplitter [22] and the same long working-distance objective lens, which is 
used for collecting the scattering radiation, or at an oblique angle using an auxiliary long 
focal-length lens. The scattered radiation collimated by the objective lens passes one or 
two notch filters (Kaiser Optical Systems) and is focused into a confocal pinhole, which is 
conjugated to the entrance slit of the spectrograph, to reduce the unwanted fluorescence 
from the diamond anvils. The best performances have been realized with the 460-mm 
focal length f/5.3 imaging spectrograph (Jobin Yvon HR460) equipped with two 1800 
and 300 grooves/mm gratings on the same turret. This feature allows one first to take 
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Fig. 1. - Layout of an optical system used for the study of metals under pressure. 



a quick spectrum over a wide spectral range (about 5000 cm^ 1 ) and then study regions 
of interest with higher accuracy and resolution. The system can be tuned to other exci- 
tation frequencies by changing the laser wavelength and rotating the holographic filters. 
One set of holographic filters covers the spectral range about 70 nm (notch filters) and 
150 nm (bandpass filter), which is suitable for most applications. Use of notch filters 
for laser radiation rejection limits the spectral range available for Raman measurements 
to Au >50-100cm _1 because the bandpass of the notch filter cannot be made narrower 
with present technology. If the low-frequency range is needed, then the spurious laser 
radiation can be supressed by a double tunable filter made with dispersive holographic 
transmission gratings [22]. 

Use of the double spatial filtering (one for the laser and one for the signal) effectively 
suppresses laser plasma lines and unwanted Raman/fluorescence signals (e.g., from dia- 
mond anvils). In this configuration the laser spot can be made very tight (e.g., 2 ^.m and 
the depth of focus is substantially reduced. Thus, the Raman signal from the diamond 
anvil is suppressed to a minimum; the effects of pressure unhomogcneity can also be 
substantially reduced [23]. Diamond fluorescence can be further reduced by using syn- 
thetic ultrapure diamonds anvils. This is crucial for studies of very weak scatters (e.g., 



4 



ALEXANDER. F. GONCHAROV ETC. 



high-temperature superconductors) and/or materials at very high pressures (above 200 
GPa), where strong stress-induced fluorescence of diamond anvils is a major obstacle for 
obtaining Raman/fluorescence spectra of samples [23]. 

In case the of metals, the use of an angular excitation geometry was found to be 
crucial, since it allows a substantial reduction in background compared to the pure back- 
scattering geometry. For this, the diamond cell must be modified to allow the angular 
light access. Specially designed tungsten carbide seats with angular conical holes have 
been used for this purpose. 

2'2. Materials. - Samples of iron-nickel alloys Fe^^Ni^ with x<0.2 from various 
sources have been studied recently. Some of them (with x=0.05, 0.10, 0.15, 0.20) were 
synthesized and characterized at Cornell University [24]. Another group with x=0.05 
and x=0.1 were the same materials studied in Ref. [25]. New samples with x=0.01, 
0.03 and 0.07 were synthesized at the Geophysical Laboratory. The sample with x=0.07 
was enriched by 57 Fe isotope for nuclear resonance inelastic x-ray scattering (NRIX) 
measurements. Rhenium samples were either commercially available polycrystallinc foils 
or powders in the case of our in hydrostatic experiments. Samples of Mg 10 B2 were similar 
to those used in Refs. [27, 28]. They are essentially in a powdered form and consist of 
aggregates of 30-50 /urn linear dimensions, which is ideal for high-pressure experiments. 

3. Raman studies of metals 

3T. Iron and Fen_ x \Ni x alloy. - Knowledge of the physical properties of iron and 
its alloys is very important for understanding the nature of the Earth's core. Recently, 
there has been a flurry of experimental and theoretical activity devoted to determination 
of the elastic properties of iron at the conditions close to those of the core [29, 30, 
31, 32, 33, 34, 35, 36]. Despite these efforts, direct and unambiguous measurements of 
the second-order elasticity tensor (e.g., CV,) under pressure are missing, and existing 
experimental and theoretical data are still controversial [31, 37]. Raman spectroscopy 
can provide independent constraints on the high-pressure elasticity from measurements 
of the frequency of the transverse optical phonon, which is related to the shear elastic 
constant C44 by the relation 

(1) C44 - 27r 2 M[V3c/(6a 2 )]^ 2 , 

where M is the atomic mass of iron, c and a are unit cell parameters, and v is the 
frequency of the optical phonon [37, 38]. This simple relation reflects the fact that the 
frequency of the phonon and the sound velocity (or corresponding Cij ) are the properties 
of the same phonon branch (fig. 2), folded because of the doubling of the fee unit cell 
to form hep. Figure 2 shows the dispersion curves in the TA direction (001) for the e-Fe 
computed theoretically [35, 39] in comparison with a "sine" curve obtained with the same 
phonon frequency. They are very close to each other (but a small discrepancy near the 
zone boundary), which show the validity of the model used. The theoretically computed 
C44 is 232 GPa, while the model gives a comparable value of 258 GPa. In general, this 
method gives fairly good results in for metals that exhibit a "normal" behavior under 
compression, while the proximity of phase transitions may cause a substantial anisotropy 
of elastic properties, which makes significantly reduced the validity of the approach (see 
discussion below). 
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Fig. 2. - Illustration of folding of the dispersion curves in a hexagonal closed-packed (hep) 
structure due to a doubling of the unit cell (with respect to a cubic one). The zone boundary 
transverse acoustic mode propagating along [001] direction becomes the zone center E2 9 Raman- 
active mode. The gray solid and dashed (folded) lines are first-principles calculations at about 
60 GPa [35, 39]; the black solid line is a "sine" curve computed with the same zone-boundary 
frequency (279.5 cm -1 ). 



We performed the Raman measurements in e-Fe to 152 GPa, which correspond to 
pressure at the core-mantle boundary [20] (see also Ref. [37] with similar measurements 
to lower pressures). Raman spectra of iron at selected pressures are shown in (fig. 3). 
The strong band, which appears in the spectra above 15 GPa on pressure increase and 
remains to approximately 10 GPa on pressure decrease, is identified as the E 2g mode 
of e-Fe. The difference for compression and decompression is consistent with the known 
hysteresis through the a-e transition [40]. A second band, which is weaker and broader 
than the main one, is observed at higher frequencies at moderate pressures. It will be 



6 



ALEXANDER F. GONCHAROV ETC. 




Fig. 3. - Raman spectra of Fe at different pressures to 152 GPa [20]. Spectra at 7 GPa and 
14 GPa are measured on pressure release. The points are experimental data; the solid lines are 
phenomenological fits (Lorentzians). The spectra are shifted vertically for clarity. 



discussed below in comparison with analogous behavior found for the alloys. 

We also present preliminary results of a similar study of Fe( 1 _ a ,)Ni a; alloys and compare 
the results with pure Fe. Ni can be considered a secondary element in the Earth's core, 
so we addressed possible change of elastic properties of the chemical substitution. The 
phase diagram of the low-nickel Fc^^Ni^ alloys are similar to pure Fe [24], and the 
transition pressures at room temperature are very close. Moreover, the equations of state 
are very similar in both the a or e phases [25, 26]. The Raman spectra of Feo.93Nio.07 at 
different pressures are shown in fig. 4. As for pure Fe, the dominant feature of the spectra 
is the E,2g mode, which appears when pressure exceeds the a—e transition boundary (>15 
GPa upstroke in agreement with numerous x-ray data). The second broader and weaker 
peak is more pronounced in the alloy compared to pure Fe (fig. 4) . It could be observed 
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Fig. 4. - Raman spectra of Feo.93Nio.07 at different pressures to 60 GPa. Spectra at 7.5 GPa 
and 15 GPa are measured on pressure release. The spectrum at 15 GPa is shown along with a 
manual fit using a coupled oscillator model [42, 43, 44]. The spectra are shifted vertically for 
clarity. 



to almost 60 GPa (c.f., 45 GPa in pure Fe) and also after the back transformation to 
a-Fe (<8 GPa downstroke), where the E2 9 mode, which is only active in e-Fe, becomes 
unobscrvable. The E2 9 mode line profile appears to be asymmetric at low pressures, 
which may be related to coupling to other excitations. 

The frequency of the E2 9 mode is lower by some 10 cm" 1 for the alloys compared 
to pure Fe, and this difference decreases with pressure, so the curves merge above 100 
GPa (fig. 5). Surprisingly, this Raman frequency offset is independent of composition 
for x>0.03. For x<0.03 the pressure versus frequency data are undistinguishable within 
the experimental accuracy. The pressure dependence of the frequency of the second 
broad Raman peak for the pure Fe and the alloys are close to each other (fig. 5). In 
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Fig. 5. - Pressure dependence of frequencies of the Raman excitations for pure Fe and repre- 
sentative Fe:Ni alloy (Feo.93Nio.07)- The curves for other x (Ni compositions) coincide within 
experimental uncertainty with either that for pure Fe (if x<0.03) or for the Feo.93Nio.07 com- 
pound shown (if x>0.03). Circles correspond to the results of fitting two-Voigt spectral functions 
(Fano plus Voigt for the 15 GPa spectrum) to the measured data. Diamonds and dotted line 
are theoretical results from Refs. [31, 45] and Refs. [35, 39], respectively. 



the case of the alloy, it seems to show an upturn with pressure. This second Raman 
peak was tentatively assigned in [20] to disorder-induced phonon scattering. According 
to calculations [41], the LA phonon branch is very flat between K and M points of the 
Brillouin zone, thus giving a large density of phonon states close in frequency to the 
experimentally observed second Raman peak. It can also be due to TA zone boundary 
phonons at the L point of the Brillouin zone [47]. The large density of phonon states in 
this frequency range is also confirmed by direct measurements [35]. The reason for the 
disorder may be a proximity to the transition boundary which can give rise to stacking 
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Table I. - Comparison of the measured and deduced from the Raman measurements shear elastic 
modulus C44 at ambient conditions. 



Element 


Raman Shift(cm x ) 


C44(GPa) calculated 


C44 (GPa) experimental 


Zr 


87 [15] 


28.9 


33.4 [50] 


Mg 


122 [17] 


15.9 


18.6 [50] 


Zn 


71.5 [18] 


19.8 


39.6 [50] 


Re 


119 (this work) 


133 


161 [51] 



faults (mixed hep- fee stacking). The alloys are also naturally disordered because of the 
presence of different atoms. In principle, disorder relaxes symmetry and wave- vector 
conservation requirements for Raman activity and causes coupling between excitations 
with different parent symmetries. The results of a coupled oscillator model fit [42, 43] 
are shown in fig. 4 [44]. An alternative explanation would be possible weak magnetism 
in e-Fe at lower pressures [31, 46], which would lift the degeneracy of the E 2s mode. Both 
interpretations are consistent with disappearance of the effect at higher pressures. 

The lower frequency of the E2 S mode in the alloys compared to pure Fe is intriguing. 
One would expect very small reduction of the frequency (0.5% or 1 cm -1 at 25 GPa 
for the Feo.sNio.2 alloy with the maximum substitution studied) if only the effect of 
changing the mass is taken into account. Since no effect of changing of lattice parameter 
with Ni substitution is reported, a substantial change in force constant is not expected. 
One possible explanation for the observed softening would be different coupling constant 
between the F,2 9 mode and the second Raman mode for the pure Fe and alloys. A larger 
coupling constant (in case of alloys with more relaxed selection rules, see above) would 
shift the E 2s mode to lower frequency. 

The threshold nature of the effect observed (see above) suggests a change in behavior 
at about x=0.03. This may be due to a change in magnetic (or chemical) ordering or a 
change of free carrier density (which in turn can influence the coupling of the F,2 g mode 
with free carriers). If qvp/hto <1, where w and q are the phonon frequency and wave 
vector, vp is the Fermi velocity, then the phonon experiences so-called Landau damping, 
so that its frequency and damping re- normalize [48]. The magnitude of the effect can 
be tuned by variation of pressure (through variation of oj) and by doping level, which 
change the Fermi velocity v^?. Additional studies need to be performed to clarify the 
observed effect. 

3'2. Rhenium. - Our initial motivation for measuring the Raman spectra of Re under 
pressure was possible use it as a convenient secondary pressure calibration gauge. As an 
hep metal, Re has a Raman-active E2 9 mode, which has been measured under pressure 
by Jephcoat and co-workers [16, 17, 18]; no anomalies have been reported consistent 
with the wide stability field of the phase. On the other hand, there are several peculiar 
observations at ambient and low pressure, which have motivated further study. Raman 
measurements at low temperatures and ambient pressure show extra peaks, which have 
been ascribed to a possible structural instability [49] . The pressure dependence of super- 
conducting transition temperature show anomalous behavior up to 1.6 GPa [53], which 
has been tentatively ascribed to a pressure-induced change of the topology of the Fermi 
surface (electronic topological transition; ETT) [54] . 

We measured the Raman spectrum of Re to 138 GPa (fig. 6). We found a positive 
monotonic shift in frequency of the E 2g mode with pressure in agreement with the results 
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Fig. 6. - Raman spectra of rhenium at selected pressures to 138 GPa. The spectrum at 138 GPa 
(dots) is shown along with Lorenzian fit (solid line). Spectra are shifted vertically for clarity. 



of Ref . [17] (which stopped at 65 GPa). A broadening of the band with pressure (fig. 6) 
arises from nonhydrostatic stresses because no pressure transmission medium was used 
(as in case of experiments reported in Ref. [17]). The frequency shift (fig. 7) is weakly 
sublinear with pressure (c.f., Ref. [17]). 

We also performed measurements under quasihydrostatic conditions to 44 GPa (using 
Ne as a medium [55]) and found that there is a systematic deviation in the frequency 
to the higher energy side compared to nonhydrostatic data [56, 57]. As the result, the 
pressure dependence of the E2 9 mode frequency under quasihydrostatic conditions is 
essentially sublinear at low pressures (fig. 7). The volume dependence of the mode- 
Griineisen parameter (inset to fig. 7) obtained from these data and x-ray diffraction 
measurements of the equation of state [58] shows a strong volume dependence in contrast 
with the measurements of Olijnyk et al. [18]. Also, the initial mode-Griineisen parameter 
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Fig. 7. - Pressure dependence of frequency of the E2 9 Raman phonon of Re. Filled circles - 
hydrostatic conditions (Ne medium), empty symbols nonhydrostatic data. The inset shows the 
pressure dependence of the mode Griineisen parameter. Solid line-our data; dotted line -Ref. 
[18]. 



7^0=2.16 is substantially larger than in Ref. [18] (7io=l-76). We relate the large 7io and 
its volume derivative to possible proximity of Re to ETT at < 1 GPa [53] (see discussion 
below) . 

Using formalism developed in Refs. [16, 20, 38] , we calculated the pressure dependence 
of shear modulus C44 for Re from our Raman measurements of the E2 9 mode frequency 
(fig. 8). The result at ambient pressure is 17% lower than the direct measurements, which 
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Fig. 8. - Pressure dependence of the shear modulus C44 of Re deduced from the E2 9 Raman 
frequency. Solid line, fit to the experimental data; squares, ultrasonic measurements [51]; solid 
circle, lattice strain measurements [52]; upward triangles, first-principles LDA calculations [31]; 
downward triangles, first- principles GGA calculations [31]. 



appears to be a systematic offset of the model based on comparisons among different 
metals (Table I). This systematic deviation can be easily explained by departure from 
the ideal "sine" shape of the dispersion curve inferred in Refs. [16, 20, 38] because of 
long-range interactions. This deviation can be substantially larger (as in the case of Zn) 
if there is a substantial departure from "normal" behavior (e.g., ETT). Nevertheless, if 
the validity of approximation [16, 20, 38] is established by an independent method at low 
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Fig. 9. - Raman spectra of the E2 9 mode in MgB2 at different pressures. Spectra are shifted 
vertically for clarity. 



pressures (by direct measurements), it can be successfully used at high pressures, which is 
very significant if other measurements and/or theoretical calculations give contradictory 
results. 

The pressure dependence of C44 for Re (fig. 8) determined in this work shows a 
systematic offset compared to first-principles theoretical calculations [31]. This deviation 
is smaller at higher pressures, and it can at least partly be accounted for by the systematic 
underestimation of C44 from the Raman frequency (see discussion below). This result 
differs from the lattice strain measurements (see discussion in Ref. [52]). 

3'3. Mgl>2 and ETT. - The recently discovered superconductor MgB2 [59] appears 
to be a clear case of a conventional superconductor with strong electron-phonon cou- 
pling [27, 60, 61, 62, 63, 64, 65, 66]. Theoretical calculations demonstrate that only one 
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Fig. 10. - Pressure dependence of the E2 9 Raman frequency [70] . The inset shows the pressure 
dependence of the linewidth. 



phonon branch corresponding to vibrations of the B atoms perpendicular to the c axis is 
responsible for most of the coupling [63, 64]. This branch becomes Raman-active at zone 
center (E 2g mode), which allows study of the electron-phonon coupling by vibrational 
spectroscopy. Application of pressure tunes the parameters of the system, so it can be 
used as an important test for theoretical models. 

A combined Raman, x-ray diffraction and T c studies under hydrostatic pressure [68, 
69, 70] to 57 GPa confirm the first-order phonon nature of the broad Raman band and 
provided the volume dependence of the electron-phonon coupling parameter. Figure 9 
shows monotonic shift of the E 2g mode to higher frequencies with pressure. This band is 
very broad because of the electron-phonon coupling [65], becoming even broader at 15-25 
GPa, which may be due to the appearance of an additional band close to the original 
one. At higher pressure the band becomes narrower, and no additional Raman band 
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Table II. - Comparison of initial mode-Griineisen parameters 7io and their volume derivatives 
qio for different metals 



Compound 


7i0 




Comment 


Reference 


diamond 


1.00(0.03) 


-1(1) 


normal behavior 


[71] 


iron (e phase) 


1.68(20) 


0.7(5) 


normal behavior 


[20] 


Zn 


3 


4.7 


close to EET 


[18] 


Re 


2.16(5) 


9(3) 


close to ETT? 


this work 


MgB 2 


2.9(3) 


0(1) 


large anharmonicity 


[68] 


MgB 2 


1.5 (P> 29 GPa) 




change of 7, is due to ETT 


[70] 



structure is observed. 

The pressure dependences of the Raman frequency and linewidth of MgB2 are pre- 
sented in fig. 10. The pressure shift of the frequency has a kink at about 20 GPa, so the 
pressure derivative of the frequency discontinuously diminishes at this point, which cor- 
responds to the step-like change in the modc-Gruneiscn parameter (Table II) because no 
anomaly is found in the P — V curve up to 50 GPa [67]. The mode-Griineisen parameter 
from ambient to 20 GPa is anomalously large, which we also ascribe to the electron- 
phonon coupling [68] and the peculiar topology of the Fermi surface [64]. It has the 
crossing of the lower B a bands with the Fermi level modulated by the zero-point motion 
associated with the E 2s mode. Pressure pushes these bands up and the amplitude of the 
zero-point motion decreases with pressure, reducing this dynamical crossing related to 
the E 2g mode. It is interesting that the T c versus pressure curve also changes slope close 
to 20 GPa [70]. 

The mode-Griineisen parameters (Table II) appear to be sensitive not only to the 
phase transformations, but also to changes in electronic structure (like ETT) through 
the electron-phonon coupling. This is due to the fact that 7i as well as T c are related 
to high-order derivative of thermodynamic potential, and thus are more sensitive to 
delicate changes in comparison to the phonon frequency or lattice parameters. Changes 
in 7i (qi) are in turn even more sensitive and may indicate the proximity to an ETT or 
an approaching instability of the lattice. 

Table II contains information about the mode-Griineisen parameters and their deriva- 
tives for several materials. Diamond is known to be stable to at least 350 GPa and can be 
considered as an example of " normal" behavior under pressure. Accordingly, 7^ is almost 
independent of pressure ji s=sl, which is typical for covalent compounds [72]. Metals are 
expected to behave similarly; this is illustrated by the example of e-Fe, which is stable to 
at least 300 GPa and exhibits no substantial change of structural parameters (e.g., c/a 
ratio) observed [26]. Zinc shows anisotropic compression [73], and various calculations 
suggest the occurrence of an ETT [74, 18], which has not been confirmed experimentally 
[18, 75] in agreement with recent calculations [76]. We speculate that the enormously 
large 7io and for Zn (and Re; fig. 7) indicate the proximity to an ETT, which could 
be at negative pressures. With its very large electron-phonon coupling, MgB 2 appears 
to be a very clear case of an ETT driven by the phonon [70] . 

4. Conclusions 

The development of optical spectroscopy instrumentation make it now possible to 
obtain the Raman spectra of metals to megabar (>100 GPa) if not multimegabar pres- 
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sures. The information obtained includes the behavior of first-order phonon modes and 
their coupling to the electronic and (possibly) magnetic excitations. It can be used 
to estimate the elastic moduli, and to obtain diagnostics of phase transformations and 
electronic topological transitions. 

* * * 
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